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In the heterogeneous cellular scenarios with macrocells, femtocells and/or picocells users may suffer 

from significant co-channel cross-tier interference. To manage this interference 3GPP has introduced 

,__i almost blank subframe (ABSF), a subframe in which the interferer tier is not allowed to transmit data. 

F Vulnerable users thus get a chance to be scheduled in the ABSFs with reduced cross-tier interference. We 

^ analyze downlink scenarios using stochastic geometry and formulate a condition for the required number 

of ABSFs based on base station placement statistics and user throughput requirements. The result is a 

^^1 semi-analytical formula that serves as a good initial estimate and offers an easy way to analyze impact of 

On network parameters on the requirements. We show that while in macro/femto scenario the residue ABSF 

psj interference can be managed well, in macro/pico scenario it effects the number of required ABSFs very 

strongly. The effect of ABSFs is subsequently demonstrated via user throughput simulations. Especially 
■^ 
^^ in the macro/pico scenario it is found that using ABSFs is advantageous for the system as the victim 



cn 



users are brought up from very low performance for the price of relatively small decrease in higher 
throughput percentiles. 



X 

C^ Index Terms 
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I. Introduction 

Almost blank subframes (ABSFs) are part of Enhanced Inter-Cell Interference Coordination 
(elCIC) framework yj that the 3GPP members have proposed [|2l as means to combat excessive 
co-channel cross-tier interference in heterogeneous network (HetNet) scenarios. HetNet scenarios 
are generally cellular network scenarios that cover different types of low-power nodes, such 
as base stations (BSs), relays or remote radio heads, as underlay to the traditional macrocell 
tier. HetNet scenarios that are specifically targeted to benefit from ABSFs are combinations 
of macrocells with closed access femtocells (macro/femto) and macrocells with open access 
picocells (macro/pico) flSl. 

Femto base station (FBS), also called Home eNodeB (HeNB), is a low-power BS that is 
deployable by the end user and connects to the core network of a cellular operator by means 
of wired broadband connection. In the eyes of the operator this is a win-win situation as the 
users (femto user equipments, FUEs) benefit from higher connection throughput while the use 
of commonly available wired broadband (such as e.g. digital subscriber line) decreases costs of 
expanding network infrastructure. In a closed access femtocell only selected users have access 
to the FBS services, creating thus a closed subscriber group (CSG). A drawback of this is that 
a non-member macro UE (MUE) that is located close to a closed access FBS can suffer from 
excessive interference. Proposed (non-ABSF) solutions for downlink interference management 
in OFDMA femtocells include for example FBS power control [4J, frequency partitioning [5|, 
O, precoding f7\, cognitive radio approach [8] and augmentation of scheduling algorithms [l9l. 

Pico base station (PBS) is practically a normal base station with lower transmit power and 
therefore smaller coverage region. The point of PBS deployment lies not in covering areas where 
macro tier signal is too low, but in augmenting the macro tier in areas where the concentration of 
MUEs is too high to be efficiently served by a macro base station (MBS). For such augmentation 
to be successful it has been shown (see [[TOll for the first suggestion) that even UEs that have 
somehow stronger signal from the closest MBS should be allowed to associate to a PBS, thus 
leading to a so-called cell range expansion (CRE) concept. In that way more UEs will be 
associated to PBSs, which leads more efficient frequency reuse and desirable traffic offloading 
from the macro tier. However, as the CRE description already suggests, some of the pico UEs 
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(PUEs) see strong interference from the macro tier. Proposed (non-ABSF) solutions for downlink 
interference management in the macro/pico scenario are for example interference cancellation 
[fTT|. frequency partitioning [12J and MBS power control [[T3l . 

The ABSF concept is based on blanking some subframes of the interferer tier and scheduling 
the especially vulnerable UEs in these subframes. The vulnerable users thus get part of radio 
resources where cross-tier interference is lower. The ABSFs are called almost blank because 
not all resource elements are allowed to be blanked - the cell- specific reference symbols (CRS) 
that are used for radio resource management (RRM) measurements and channel estimation have 
to remain present. The strong interference in CRS resource elements is a separate issue and 
was suggested to be tackled by interference cancellation (see ||T41 or [[T5l also for other control 
channel challenges), but such considerations are out of scope of this work. Alternatively to 
ABSF, the BS can configure an empty MBSFN (Multicast-Broadcast Single Frequency Network) 
subframe, but its use is more constrained, therefore, our focus will be on the ABSF. Compared 
to other mentioned interference management solutions, the ABSF concept is simple enough to 
be incorporated into often technically entangled 3GPP specifications and at the same time it has 
found rather wide acceptance among the standardization partners. 

The interferer tier in the macro/femto scenario is the femto tier, while in the macro/pico 
scenario it is the macro tier. In case there are both FBSs and PBSs within MBS coverage, 
there might be need for ABSFs in femto tier as well as in macro tier. As the (significant) 
cross-tier interference can come from multiple BSs, the amount and position of ABSFs has to 
be coordinated within the network. Indeed, the organization of ABSFs is planned as a part of 
self organizing network (SON) concept f[6\. In our work we will propose how the number 
of ABSFs for downlink interference management can be set globally based on BS placement 
statistics. Such relationship can then serve as an initial estimate or as a backup solution when 
the distributed coordination does not serve its purpose. We will derive the necessary number of 
ABSFs for macro/femto and macro/pico scenarios separately and, if needed, the results can be 
easily combined. 

To our best knowledge the question of deriving the number of necessary ABSFs has not been 
addressed in any work so far. Besides the mentioned concept research (Si, [fTSl the work is 
quite sparse. For the macro/femto scenario, some simulation results of using ABSFs have been 
published in 0171 and [fTSl . In [fT9ll the authors introduce a coordination framework for ABSFs, 
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including channel quality indication (CQI) processing, and suggest control messages that are 
needed for such operation. Simulated performance of the macro/pico scenario has been shown 
in [l20l . while [|2T| presents also analytical insight into the topic. 

We address the problem by setting the number of ABSFs globally using tools from stochastic 
geometry ll22l . [|23l . Base stations and users are modeled as 2D stochastic processes and spatial 
relation between a user and its closest interferer is leveraged to define victim users, i.e., users that 
require interference management. Parametrization of the stochastic models and properties of the 
victim users are then used to formulate the number of necessary ABSFs. We thus give a semi- 
closed form connection between the stochastic intensity and other parameters and the minimum 
number of subframes that can be quickly used to determine the fraction of radio resources 
needed for interference management. Subsequently, we analyze dependence on individual model 
parameters, the most important result from which is that in macro/pico scenario the residue 
ABSF interference has a strong effect on the required number of ABSFs. Finally, we demonstrate 
the effect of derived number of ABSFs on user throughput via simulations. The results show 
moderate performance gain for victim users in macro/femto scenario, but in macro/pico scenario 
the improvement is substantial. 

The remainder of the paper is structured as follows: Section In] introduces the stochastic 



geometry -based system model and defines victim UEs. In Section III we derive the success 



probability of victim UEs, i.e., probability that signal-to-interference ratio (SIR) is higher than 



a predefined threshold. In Section IV we use the success probability to set condition for the 



necessary number of ABSFs. In Section IV] we evaluate the effect of ABSFs by means of Monte 



Carlo simulations. Section VI concludes the paper. Most of the sections are divided into two 



parts, one for macro/femto scenario and one for macro/pico scenario. 

II. System model 

In this work we model BS and UE placements as homogeneous Poisson point processes (PPP). 
In [|24ll it has been shown that the random BS placement produces a good lower bound for 
SINR distribution, a regular grid BS deployment gives an upper bound, and the actual truth lies 
somewhere in between. In [|25ll the authors have taken two samples of real world BS placements 
and shown that PPP is not actually a good model for them, because it lacks interaction between 
points. Nevertheless, the use of PPP model is prevalent (see e.g. [|26ll for heterogeneous networks) 

April 15, 2013 DRAFT 



as it offers a rare analytical insight into a larger scale network. 

Radio channel conditions are modeled by a combination of distance dependent path loss 
H{r) = r^", where a is path loss exponent, and fast Rayleigh fading with exponential power 
distribution h ~ exp(l). For the sake of tractability we consider two general simplifications. 
Firstly, we do not model shadow fading. Although it is possible to include shadow fading in 
initial model equations, a considerable degree of tractability is lost (see [|24l and [|27l '). Secondly, 
the residue ABSF interference is considered white. In reality, the victim receiver would see full 
interference on resource elements where the interferer transmits CRS, and some leakage on other 
resource elements. However, considering that the CRS position varies between cells and all data 
is subject to scrambling, the white interference model is not extremely far-fetched. 

We note here that although 3GPP specifications and state of the art research work offer quite 
a few techniques that could supplement the use of ABSF for cross-tier interference management, 
the scope of this study and structure of the system model does not allow us to consider them. 
We assume a single antenna transceiver at both BS and UE, hence multi-antenna techniques are 
not considered. Coordinated multi-point transmission and advanced receiver processing are also 
out of our scope. 

A. Macro/femto scenario 

In the first scenario we have an overlay macrocell PPP 0mbs of intensity Am and an underlay 
femtocell PPP ©fbs of intensity Ap. Macro UEs (MUEs) form another PPP 6mue of intensity 
Amue- Because we are interested in protecting MUEs from CSG femtocells, we do not need to 
model femto UEs. Processes ©mbs, ©fbs and 6mue are all separate and independent among each 
other. Indoor/outdoor positions and walls are not considered in the model. MUE is associated to 
a macro BS (MBS) with the best long-term channel conditions, i.e., the geographically closest 
one, as in [|24ll . Transmission powers of MBS and FBS are denoted by Pu and Pp, respectively. 
Path loss exponent on MBS-UE links is «m, on FBS-UE links it is «f- Data traffic in macro and 
femto layer is modeled by base station load values 0m and ^p^ respectively. A comprehensive 
analysis of load impact on PPP-based model of cellular network has been done in J28l|. Although 
our work uses a different UE-BS association model, the traffic model philosophy remains the 
same. 

Let us now have an MUE that is associated to MBS of distance tm- We say that an FBS 
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is a dominant interferer (DI) to given MUE if the MUE-FBS distance rp is smaller than kru, 
where A; is a Dl-defining coefficient. We then define a victim MUE as an MUE that has one 
or more DIs. MUE can measure long term channel conditions (and thus estimate distance) of 
FBSs by performing a reference- signal-received-power (RSRP) measurement. Presence of DIs 
(or victim status) is then reported to MBS, which uses the knowledge to decide which MUEs 
will be scheduled in ABSF and which in normal subframe (NSF). 

Value of the Dl-defining coefficient k is important here. Its basic purpose is to take into 
account difference between MBS and FBS transmission powers and path loss exponents and 
consider the maximum interferer power that the MUE receiver can withstand. If, for example, 
«M = ap = a and the minimum required SIR at MUE is OdB, k = (Pp/Pm) ■ With om 7^ «f 
the relation between received powers can no longer be transformed to linear relation between 
rp and tm and our condition becomes only approximate. 

The definition of victim MUE gives k also another dimension. It might happen that interference 
from multiple FBSs that are not marked as DIs is unbearable and the given MUE should be 
marked as victim. A safety margin in k value is an option how to counterweight this issue. 

Our definition of DIs via k resembles [|29ll where a contour of equal power is used to decide 
how to divide available spectrum. The authors however do not pursue its effect on SIR, nor do 
they analyze it from statistical point of view. 

B. Macro/pico scenario 

In the second scenario there is again an overlay macrocell PPP 6mbs of intensity Am and 
an underlay picocell PPP ©pbs of intensity Ap. Macro UEs (MUEs) form again a PPP 6mue 
of intensity Amue- Processes ©mbs? ©pes and 6mue are separate and independent among each 
other. As we now have an open access policy, part of the MUEs will actually be associated to 
PBSs and thus called PUEs. Drawing inspiration from [26], an MUE is associated to a PBS if 
the distance to the closest PBS rp is smaller than kiryi, where tm is the distance to the closest 
MBS. In case there is no PBS within kiVyi distance, the MUE is associated to the closest MBS. 
The ki coefficient takes into account differences in MBS and PBS transmission powers Pm 
and Pp, respectively, and the association bias k. Unlike OOl that tries to find a framework for 
K optimization, we keep the value constant and analyze its effect on the required number of 
ABSFs. Path loss exponent on MBS-UE links is am, on PBS-UE links it is ap. The MBS and 
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PBS load values are denoted as 0m and 0p, respectively. 

The dominant interferers to a PUE are MBSs that fulfill rp > k2rM, where rp is distance 
between PUE and the associated PBS, k2 is the Dl-defining coefficient and ru is distance to 
given MBS. With positive association bias we have ki > ^2. Victim PUEs are PUEs that have 
one or more DIs. They are thus identified by a pair of inequalities A;2'"m <rp< kiVu- Presence 
of DIs (or victim status) are reported from PUE to PBS, which can then use the knowledge to 
request ABSFs from the macro tier. The k2 coefficient has the same meaning as k in macro/femto 
scenario and can as well be used to address the issue of multiple DIs. 

III. Success probability 
In this section we derive success probabilities of victim MUEs and PUEs in our scenarios. 



which we then use in Section IV to set the necessary number of ABSFs. We use the name success 
probability as in [|26ll . i.e., a probability that UE has a signal-to-interference-ratio (SIR) higher 
than an outage threshold. SIR is used in this work to approximate signal-to-interference-plus- 
noise-ratio (SINR), as we are modeling interference limited networks. The success probability 
represents a CCDF of SIR. While for arbitrary located UEs equivalent result have been presented 
e.g. in [|3T1l , [|32l and [|33l , conditioning on presence of dominant interferers has to our best 
knowledge not been done before. 

A. Macro/femto scenario 

Signal-to-interference-ratio at a victim MUE on single resource block is given by formula 

j= ^'^^^M ^ (J) 

-^M + Pa (-^di + -^f) ' 
where h denotes fast fading power gain, Jm denotes sum interference power from the macro 
tier (all MBSs except associated one), Jdi denotes sum interference power from DIs, Ip denotes 
sum interference power from non-DI FBSs and p^ denotes residue ABSF interference. Results 
in this section are derived for ABSF, for NSF one would simply omit pa- Success probability 
is defined as 

P{7>7o}, (2) 

where 70 represents the outage threshold. 
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To increase clarity and give insight on how the work has progressed we first derive success 
probability in case of full load and single DI present. After that, we generalize it for arbitrary 
load values and one or more DIs present. 

1) A single dominant interferer: Probability that the number of FBS DIs N^^ around a 
randomly chosen MUE equals one is given from the definition of PPP as 



P|M^^ 






DI 



1 



^M 



vr/c^Apr^exp (— Tr/c^Apr^) • 



By averaging over tm we get 
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By differentiation we get a PDF 

fru\kir) = 27r' (A:'Ap + Am)' r' exp {-n {k^X^ + Am) r') . 



(3) 



(4) 
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The success probability of a victim MUE is 



P 7>7o 
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r(F) 



The reader may notice that on RHS of ( [T3| ) we have omitted the condition A^^j'^ = 1. This is 
purely for space purposes and we will repeat it a few times throughout the paper. The inner 



probability from ( [13] ) can be found as 
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The first and the last terms in (fTTl) have been derived in Il24l and are given by 



E/m {} = ^^P (-TrAM^V (7o, «m)) 



and 



where 



rWn = Pvnf-^^2,..2,f70PA^R^ 



E'^^^ {} = exp -TTk'Xpr'p 
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1 
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(14) 
(15) 
(16) 
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For averaging the dominant interference term in ( [T7| ) we exploit the Laplace transform of 
exponential function as given here 

ES{exp(-./Di)} 



E 



S {exP (-s^F/iDiror) } 

1 



e(^) , 



kr 



1 + sPpU-"'^ 

1 
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(21) 
(22) 

(23) 
(24) 
(25) 



where /iqi and tdi are are fast fading power and link distance between victim MUE and DI FBS, 



respectively, and 2F1 () is the hypergeometric function. For s- 



Pu 



this gives 



E^ < exp 
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The PDF of tm ( [12] ) and interference terms ( [T8| ), ( [T9| ) and ([26]) can now be plugged into ( [T3[ ). 
In case of single path loss exponent aM = «F = « the integration variable disappears from inside 
the p{) and 2-^1 () functions and we get a closed form solution 
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2) One or more dominant interferers: A probability that an MUE has one or more FBS 
dominant interferers within kry[ distance is complementary to probability that there are no FBSs 
within that distance. Using the same approach as with the single DI we get the probability 



P <! N^^ > 1 



A;2Ap + Am 



(28) 
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and the PDF of distance from closest MBS 
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The success probability is again given by 
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The Jm and Jp terms generalized for arbitrary BS load values are given here: 

JE/m {} = exp (-tt^mAm^V (7o, "m)) 
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The Jdi term from the single DI case can be augmented based on the independence of FBSs, 
resulting in an approximation 
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where A^^j is the average number of PBS DIs within kru distance biased by knowledge that 



— (F) 

there is at least one. The value of A^^j is 



A^- 



(F) 



vr/c^Apr^ 



(35) 



""' 1 - exp i-nk^Xprl,) ' 
where the numerator represents the unbiased average and the denominator represents the bias, 
i.e., removing the probability that there are no DIs. 
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The interference terms can now be put into pO]) to calculate the success probability. Unlike 



p7] ), the integral in ( [30| ) cannot be simplified into a more digestible form even with single path 
loss exponent and has to be evaluated numerically. 



B. Macro/pico scenario 

In the second scenario, the downlink signal-to-interference-ratio at a victim PUE is defined as 

ip + Pa (jdi + Jm) 
where Jp is sum interference power from the pico tier (all PBSs except associated one), Jdi is 
sum interference power from dominant MBS interferers and Jm is sum interference power from 
all other MBS interferers. In the rest of the subsection we will state the most important results 
for the success probability. Derivation follows the same logic as in the first scenario. 
The probability that a UE is actually a victim PUE is 

Am Am 
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and the PDF of rp of a victim PUE 
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The separate interference terms are given here: 
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where A^^j is the average number of MBS DIs calculated as 
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Like in the first scenario, even with single path loss exponent the final integral (|39]) can be 
evaluated only numerically. 

IV. Number OF ABSFs 

3GPP is introducing ABSFs to protect victim UEs. In both considered scenarios, victim UEs 
as we defined them are easily identified by presence of dominant interferers in the vicinity based 
on RSRP measurements. We can thus focus all our attention on the victim UEs and avoid a 
general and comphcated sum rate optimization problem with a simphfied one 

maximize Cwv (Na) 



subject to Cv (Na) > Cv.min, 
< A^A < A^s, 



(45) 



where A's represents number of subframes in a frame, A^a represents number of ABSFs in a 
frame, Cnv stands for throughput of non-victim UEs, Cy stands for victim UE throughput and 
Cvmin stands for the minimum required victim UE throughput, a parameter of choice. We will not 
define Cnv (^a) more closely because any reasonable definition is a strictly decreasing function 
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of JVa- We thus turn even more attention to the victim UEs and define a greatly simplified 
problem 

minimize Na 

subject to Cv (Na) > Cymin, (46) 

0<Na< Ns. 

Now, the throughput of victim UEs depends on quite many things, from which the most important 
one is how the BS schedules the UEs, i.e., what is the scheduling algorithm and how are victim 
and non- victim UE transmissions placed in NSFs and ABSFs. In [ITTll it has been suggested that 
proportionally fair scheduler, assuming reliable knowledge of channel state, can take care of the 
division between NSFs and ABSFs itself. However, scheduling is (and will probably long stay) 
implementation specific, so that BS hardware and software vendors can compete with each other 
and fight for the technological edge. 

Therefore, to define Cy (^Va), we consider what is from the throughput and scheduling point 
of view the worst case scenario: average throughput at outage threshold (outage throughput) and 
round robin scheduling algorithm. Concerning NSF and ABSF division we assume that non- 
victim UEs have access only to NSFs while victim UEs have access to both NSFs and ABSFs. 
With these the average victim UE outage throughput is 

Cv = E,,,. {^Ci^)(L,Lv) + ^^^Cr(L,Lv)} (47) 

where Cy is average victim UE outage throughput in ABSFs, Cy is average victim UE outage 
throughput in NSFs and L and Ly represent the number of all UEs and the number of victim 
UEs associated with the BS, respectively. Although L and Ly are correlated, for our purposes 
the outage throughput Cy can be very well approximated by considering them as independent. 
We thus write 

C,^^Ci'\L,) + '^^^C?\L). (48) 

iVs iVs 

Now, assuming channel independence across resource blocks and using the success probabilities 



derived in Section 



III 



we can approximate Cy and Cy as 

Ci^\Ly) ^ ATrP {7(^) > 7o} log (1 + 7o) ^^^\ (49) 

^^(L) ^ atrP {7W > 70} log (1 + 70) n^''\ (50) 
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Figure 1. Approximate outage throughput of victim MUEs in macro/femto scenario versus DI definition via e (as in k 



[Pp/ (ePm)] 



2/(aM+«F) 



) compared to simulated equivalent. 
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Figure 2. Dependence of required number of ABSFs A^a in macro/femto scenario on FBS intensity Af with other parameters 
constant. The Af/Am ratio equals mean number of FBSs per MBS coverage. 



where A^r is number of resource blocks, the fi and fi terms denote the average asymptotic 
proportion of resources that a victim UE is scheduled via the round robin principle in ABSF 
and NSF, respectively, and 7^^-' and 7*^^^^ denote SIR in ABSF and NSF, respectively. In place of 
the success probabilities we use the corresponding macro/femto and macro/pico derived terms 



t(a) 



7(N) 



( [30| ) and ([39]), respectively. The Q and Q values depend on the average number of victim 
and non-victim UEs per BS. Derivation of these values is presented in the Appendix. Finally, 
we connect all the acquired results and construct a condition for the number of ABSFs as 



N. 



mm 



-^S ( C'vmin — Cy 



^(A) 






(N) 



A^s 



(51) 



A. Macro/femto scenario 

Before moving on to performance evaluation, we demonstrate precision of our outage through- 



put approximation and analyze dependence of the condition ( [51] ) on key parameters. Unless 
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Table I 
Reference simulation parameters 



Parameter 




Value 


Simulated area 


10000m X 10000m 


Samples collected from 


3000m X 3000m 


Macro BS intensity Am 


lO-^m-2 


Femto BS intensity Af 




12Am 


Pico BS intensity Ap 




4Am 


Macro UE intensity Amue 




20Am 


Macro BS transmission power Pm 




43dBm 


Femto BS transmission power Pf 




20dBm 


Pico BS transmission power Pp 




30dBm 


Macro BS load (jtu 




1 


Femto BS load 4>f 




0.5 


Pico BS load 0? 




0.8 


MBS-UE path loss exponent qm 




2.5 


FBS-UE path loss exponent qf 




3.5 


PBS-UE path loss exponent Qp 


(ft) 


3 

2 


Macro/femto Dl-defining k 


"M+=F ^0.136 


Macro/pico association bias k 


(ft) 


7dB 

2 


Macro/pico association-defining fci 


2 


Macro/pico Dl-defining k'z 


=M+=P ^0.262 


ABSF residue interference pa 




-20dB 


Outage threshold 70 




-5dB 


Number of subframes A^s 




10 


Number of resource blocks A^'r 




25 


Resource block bandwidth 




180kHz 



Stated otherwise, parameters are at their reference values as shown in Table |Ij We consider 
these to be realistic values. Because aMT^ap and auT^ctp, we decided to use (ttM+^p) /2 and 
(ttM + ttp) /2 in the definitions of k, ki and k2. That way we keep the definitions close to the 
intuitive shape with single a. Impact of different Dl-defining coefficients is later considered in 
the performance evaluation. Unless stated otherwise, the minimum average outage throughput 
Cv.min in macro/femto scenario is set to 40kbits/s. 

In Fig. [T] we plot the approximation of outage throughput compared to simulated equivalent. 
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Figure 3. Dependence of the required number of ABSFs Na in macro/femto scenario on residue ABSF interference pA in 
macro/femto scenario with other parameters constant. 



Macrocell path loss and minimum average outage throughput are changed from default values 
to au = 3 and Cymin = 60kbits/s in order to increase simulation precision (for further discussion 
see Section [V]). Other parameters are at their default values as in Table |lj The values of Cy are 
plotted against definition of DI. The horizontal axis does not present k directly, but a ratio e that 
puts into relation own received power and DI received power, i.e., k= [Pp/ {^Pu}] '^°'^~^^'. The 
approximation can be considered very good, with deviations coming mostly from insufficient 
number of Monte Carlo samples. Within e E (— 4dB, — 3dB) and e G (2dB, 3dB) there is a step 
in the curve because for given values of k the number of required ABSFs according to our 



condition (51) changes. 



Dependence of ( [51] ) on Cymin and Amue is intuitively clear and our results confirm it - the A^a 
required grows approximately linearly with both of these values. Therefore, we save space and 
exclude those results. In Fig. |2] we show dependence on Ap. The trend is approximately linear 
in the beginning and then becomes slightly saturated. In Fig. [3] we show dependence on pa, 
the residual FBS interference in ABSF. Here the dependence is low because FBS interference 
is strongly attenuated by larger path loss exponent ap. Finally in Fig. |4] we show dependence 
of Aa required on the definition of dominant FBS interferer k via e. With increasing e the k 
coefficient is decreasing and less MUEs are considered being victim. Although this decreases 
the Aa required, we are practically increasing the number of non-victim MUEs and therefore 
have to consider effect on their performance. We will come back to this in Section |Vl 
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Figure 4. Dependence of the required number of ABSFs A'a in macro/femto scenario on the definition of victim MUE with 
other parameters constant. 




2 4 6 8 10 12 

PBS intensity/MBS intensity 



Figure 5. Dependence of the required number of ABSFs A'a in macro/pico scenario on PBS intensity Ap with other parameters 
constant. The Ap/Am ratio equals mean number of PBSs per MBS coverage. 




Figure 6. Dependence of the required number of ABSFs A'a in macro/pico scenario on residue ABSF interference pA with 
other parameters constant. 



B. Macro/pico scenario 

Similarly to the previous subsection, we present dependence of A^a on selected parameters. The 
reference parameter values are taken from Table |I] and the minimum average outage throughput 
Cvmin is set to lOOkbits/s. 

In Fig. [5] we show the effect of Ap value on A^a with other parameters at reference values. 
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Figure 7. Dependence of the required number of ABSFs A'a in macro/pico scenario on association bias k, which defines fci, 
with other parameters constant. 



Compared to intensity of FBSs in macro/femto scenario, increasing number of PBSs does not 
primarily increase interference. On the contrary, with higher number of PBSs the MUEs have 
more BSs to connect to, therefore number of PUEs per PBS decreases and less ABSFs is needed. 
In Fig. [6] we show dependence of Nj^ on residue ABSF interference Pa- While below Pa = — 15dB 
the residue interference seems manageable, higher values leads to dramatic increase in the A^a 
requirement. The impact is much stronger than in macro/femto scenario. It is partly because of 
lower path loss exponent on MBS-UE links au < ap and partly because of the vulnerability 
imposed by the association bias k. Finally in Fig. |7] we present dependence of Aa on k. As 
intuitively expected, increasing k results in MBS interferers being even closer to the victim PUE 
and thus the Aa requirement increases. 

V. Performance analysis 

In this section we demonstrate the effect of the derived number of ABSFs on UE throughput 
in downlink. We use Monte Carlo simulations to evaluate UE throughput with link adaptation 
modeled by Shannon's capacity. The simulations consist of snapshots, during which BSs and 
UEs are dropped and kept at fixed positions, and in each snapshot there are multiple frames 
(with subframes) where the UEs are scheduled. Reference simulation parameters are summarized 
in Table |lj Although we consider the model parameters realistic, the computational complexity 
of the simulations prevented us from using a simulated area size that would provide sufficient 
precision for au = 2.5. However, rather than using unrealistic parameters we keep the area 
at manageable value and accept that the performance results in this section are a little on the 
optimistic side. As a measure of performance we collect aggregate throughput of each MUE/PUE 
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Figure 8. Throughput CDF of victim MUEs in macro/femto scenario with round robin and proportional fair scheduling, with 
and without ABSFs. The required outage throughput is marked by a dashed line. 
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Figure 9. Throughput CDF of all MUEs in macro/femto scenario with round robin and proportional fair scheduling, with and 
without ABSFs. 



and then evaluate the resuks in form of UE throughput CDF. 

A. Macro/femto scenario 

The minimum average outage throughput in macro/femto scenario is set to 40kbits/s. In Fig. 
[8] we present throughput CDFs of victim MUEs. Victim MUEs are scheduled in both NSFs and 
ABSFs and for comparison we have included both round robin and proportionally fair scheduling 
algorithms. The proportionally fair scheduler requires channel knowledge and in each RB it 
chooses the user with highest instantaneous SIR normalized by its mean SIR. The dashed line 
in Fig. [8] represents the chosen minimum mean victim outage throughput Cymin- 
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Figure 10. Mean round robin tlirougiiput of victim/non-victim/all MUEs in macro/femto scenario versus definition of dominant 
interferes 

Table 11 
ViCTIlVI MUE OUTAGE THROUGHPUT AND NON-VICTIM MUE IVIEAN THROUGHPUT IN KBPS FOR DIFFERENT NUMBERS OF 

ABSFS IN MACRO/FEMTO SCENARIO. 



Na 


1 


2 


3 


4 


5 


Victim mean outage TP 


32 


43.6 


53.8 


65.3 


75.4 


Non-victim mean TP 


188.2 


170.7 


148.7 


128.1 


107.7 



Looking at the victim MUE curves we can observe that in given macro/femto scenario the 
effect of ABSFs is rather moderate. Without them, around 70% of victim MUEs have lower 
throughput than Cymin- With ABSFs, approximately 50% have higher throughput than Cymin- 
With proportionally fair scheduler the throughput values are higher and effect of ABSFs weaker. 
The large difference between our round robin based requirement and the proportionally fair 
performance suggests that our research should be expanded by considering advanced scheduling 
during the ABSF planning phase. The effect of ABSFs on purely non- victim MUEs is not shown, 
as those obviously lose 30% of available resources. In Fig. |9] we show throughput CDFs of all 
MUEs combined. The rather small improvement of victim MUE performance seems to be in 
overall statistics completely overshadowed by the effect of resource restrictions for non-victim 
users. From the system perspective these findings suggest that rather than blocking femto layer 
users and non-victim MUEs from a fraction of resources, the MUEs that suffer from strong FBS 
interference should be treated in an individual manner. 



In Fig. 10 we present impact of the Dl-defining coefficient k to MUE throughput via the same 
e factor as in Subsection HV-A[ Because differences in CDF curves are too small to notice with 
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Figure 1 1 . Throughput CDF of victim PUEs in macro/pico scenario with round robin and proportional fair scheduling, with 
and without ABSFs. The minimum average outage throughput is marked by a dashed line. 



eyes, we put measured mean throughput on the y-axis. In an interval with the same number of 
ABSFs increasing value of e leads to less victim and more non-victim MUEs per MBS, therefore 
victim MUE mean throughput is increasing, while non-victim MUE throughput is decreasing. 
For our parameters a good point of operation is with e E (2dB,3dB), where the number of 
ABSFs is lower than around £ = — 2dB and the mean victim and non- victim MUE throughputs 
are relatively close to each other. 

To justify the simplification of optimization problem ( |45| ) into ( |46l ) we present in Table [ll] 
mean outage throughput of victim MUEs and mean throughput of non-victim MUEs for different 
numbers of ABSFs A^a- Mean throughput of non-victim MUEs is a clearly decreasing function 
of ABSFs, the simplification is thus well justified. Although the table shows that 2 ABSFs would 
suffice to fulfill Cymin = 40kbps requirement, this is because our simulation area is not large 
enough as mentioned at the beginning of this section. 



B. Macro/pico scenario 

In the second scenario we set the Cymin value to lOOkbits/s. In Fig.[TT]we present throughput 
CDFs of victim PUEs with the same round robin and proportionally fair scheduling algorithms. 
As before, victim PUEs are scheduled in both ABSFs and NSFs, while non-victim PUEs have 
access only to NSFs. Without ABSFs, the victim PUEs experience very bad performance. 
Strong cross-layer interference (with au < ftp) with association bias k on the top leads to 
very low throughput values. With both round robin and proportionally fair scheduler, more than 
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Figure 12. Throughput CDF of all PUEs in macro/pico scenario with round robin and proportional fair scheduling, with and 
without ABSFs. 



95% of victim PUEs experience lower throughput than the Cymin requirement. With ABSFs, 
approximately 70% and 90% of victim PUEs reach the requirements using round robin and 
proportionally fair scheduler, respectively. 



In Fig. 12 we plot results from all PUEs. Because of the association bias the proportion of 
victim PUEs is quite high, therefore the decrease of non-victim PUE performance is not even 
visible in the figure. Overall, even though there are still PUEs experiencing zero performance. 



ABSFs substantially improve service in the pico layer. For completeness. Fig. 13 also shows 
throughput CDFs for all UEs in the scenario, i.e., PUEs as well as MUEs. It confirms that ABSFs 
are very advantageous for given scenario - high performance gains in the lower percentiles are 
balanced by relatively small performance decrease in the higher end of the curve. 



A sharp reader looking at Fig. 12 and Fig. 13 notices that despite using ABSFs, a percentage 



of PUEs is experiencing system outage, i.e., no throughput at all. These are PUEs that did not 
fulfill the DI definition via reference k2. To further investigate this important design property 



we plot in Fig. 14 throughput CDFs of all UEs in macro/pico scenario with multiple values 
of e that defines k2 = [Pp/ (ePu)] '^"^^^^K Decreasing the DI definition threshold noticeably 
increases throughput values in the lower percentiles and thus reduces system outage. Because 
our reference scenario settings yield relatively low amount of PUEs per PBS, increasing the 
number of victim PUEs actually does not show any penalty in this result set. This will however 
not hold in general, as the PBSs are expected to be operational especially in areas with higher 
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Figure 13. Throughput CDF of all MUEs and PUEs in macro/pico scenario with round robin and proportional fair scheduling, 
with and without ABSFs. 
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Figure 14. Throughput CDF of all MUEs and PUEs in macro/pico scenario with ABSFs, round robin scheduling and several 
values of Dl-defining k2 via e. 



UE density. 

VI. Conclusions 

Almost blank subframes (ABSFs) offer a simple and efficient way of decreasing the level of 
background cross-tier interference and thus give an opportunity to serve vulnerable users. We 
propose a way to approximate the required number of ABSFs based on Poisson point process 
network deployment statistics. We derive the necessary number of ABSFs as a formula that is 
easy to evaluate for macro/femto scenario with closed subscriber groups and macro/pico scenario 
with cell range expansion. We analyze dependence of the result on individual parameters, showing 
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^AMUE-16807AMUE(y^^;;^)' 






7^(A) 



'"" ^ ^m\ 343SI/' (754 + 2Amue)'/' ^ (75* + 4Amue)'/' / V 3435^^" (75* + 2Amue)'/' , 



that while in macro/femto scenario the white residue interference in ABSF can be tolerated well, 
in macro/pico scenario its effect on the required number of ABSFs is substantial. Throughput 
simulations show that in macro/femto scenario the performance gain when using ABSFs is 
rather moderate and the system may work better with a more individual treatment of high FBS 
interference. On the other hand, in the macro/pico scenario we can see that, especially due to 
the association bias, using ABSFs improves the performance of the system considerably. 

The fraction of victim users in outage decreases from over 95% to 30% or even 10%, depending 
on the scheduling algorithm. Looking at all users in the macro/pico scenario, the small decrease 
in high throughput percentiles is more than balanced by substantial performance increase in the 
lower region. 

Appendix 
Derivation of round robin resource fractions 

If a loaded BS governs L associated UEs and all UEs have uniform traffic requirements, 
each UE will asymptotically be scheduled in 1/L fraction of available resource blocks. We are 
interested in expected value of this fraction for victim UEs. Using a clever approach from [;34i 
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we can write 



— °° 1 r 1 1 

^ = CnY.fW = -f\ (55) 

oo ^ 

= Cn5^/P{L = £} (56) 

OO „QQ 

= CnJ2 ^{L = i\A = x}fAix)dx (57) 

oo „ 



^^^e-^--/^(x)da; (58) 



= C^ /°° e-^--A(x) 5^ ^^^^dx (59) 

= Cn / (1 - e-^""'^^) fAix)dx, (60) 

where /^(a;) is a PDF of an area where the given set of UEs can be located and C'q and 
Cq are normahzation constants. We denote the area PDF /a(x) as fljf (x) for all MUEs in 
macro/femto scenario (round robin fraction in NSF), f^f{x) for victim MUEs in macro/femto 
scenario (round robin fraction in ABSF), f^„{x) for all PUEs in macro/pico scenario and f^j^ix) 
for victim PUEs in macro/pico scenario. The normalization constant C'q is needed because we 
are excluding cases with no UEs associated to BS, while Cq takes also into account multiplying 
the PMF F{L = i} with a weight factor i and can be calculated from 

oo 

C^' = Y,iF{L = i} (61) 



e=i 



/•oo 
/ XMVEXfA{x)dx. (62) 



'0 

In macro/femto scenario the PDF f^^^ (x) can be obtained from the general approximation of a 
Poisson Voronoi cell area that has been found in [l35l (and used in [[34l|) by adding a condition 
that there is at least one victim MUE present. The general Voronoi area for our macro PPP is 
drawn from PDF 

fs{x) ^ Am — J— (AMa;)2exp (--Xuxj . (63) 



April 15, 2013 DRAFT 



27 



The PDF f^^' (x) can be obtained from CDF 

Fs{x,p^'^^)=F{S<x\Ly>0} (64) 

_F{S<x,Ly>0} ^^^^ 

(66) 
(67) 

where p^^^ is a probability that a UE is a victim UE, i.e., in macro/femto scenario p^l = 
P \ NJ^i ^ 1 f from ( [28] ). After calculating the integrals in ([68]) and differentiating the CDF we 



F{L 

/;p{i 


v>0} 
.v>0|5 = 


= u} fs{u)du 


/;(i-p{Lv= 


= u} fs{x)dx 
0\S = u})fs{u)du 




F{Ly = 0\S = u})fs{x)dx 



get a result ( [52) ) with subterms: 

So = 7Am + 2Amue (69) 

r, _ k AfAmUE ,r.r.~. 

fc^Ap + Am 
^2 = (7A^ + A;%So) (71) 



_ /Am(/c2Af + Am) ,„, 

^^ - V 7A^ + PAp^o ^^^^ 

The PDF flajf{x) of area where victim MUEs can be located is not trivial to find. The most 
sensible approach we came up with is transforming the unconditional Voronoi cell random 
variable using again the probability that MUE is a victim MUE, i.e., 

Pmlf \Pmlf / 

In macro/pico scenario the pdf of area of all UEs within a macrocell conditioned on presence 
of a victim PUE is 

Ui^) = j^Fs{xJ^l), (74) 



where p)^ =¥ {k2rM < t? < kiryi] as in ( |37] ). For PUEs and victim PUEs we make a similar 



approximation as in (73 ) with addition of further multiplying the transformation coefficients with 
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the average number of PBSs within MBS coverage, resulting in 

(Af) / X _ -^M 1 



/m/p {^J - ApP{rp<fcirM} 
„ , Am X 

^ Jm/p 



ApP{rp<A;irM} 
where 



(75) 



^lAp + Am 
and 



P{rp<A;irM} = 7^T^^^ (76) 



/m/p (^) ~ \ (V) •^ "^P \ (V) ■ ^ -* 

^PPm/p V^PPm/p/ 

The round robin fraction values Q can be obtained in closed form for all cases. For macro/femto 



scenario we present them in (53) and (54) with subterms given here: 



Formulas for fj^p{x), fj^/pi^)^ ^m/p ^^'^ ^m/p we leave out of the paper because they are very 
spacious and do not have enough added value by themselves. To illustrate the precision of our 



approximations we show the il values as functions of Amue/Am in Fig. 15 The approximation 



— (N) 

of Q^f is very good. In other cases our derived formulas give consistently lower values than the 
simulations. The reason behind this is that when transforming the Voronoi area PDFs fm/fix) and 
/m/p(x) we neglect that dangerous zones around FBSs and coverage areas of PBSs can overlap. 
Nevertheless, the derived values fi^^f, i^^^p and fi^/p serve as a good lower bound and ensure 



that the required number of ABSFs as given in Section IV will not be too low. 
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